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1. Introduction 



The search for ways to solve problems in the quantum field theory with a large coupling 
constant, has long been significant. Calculations on space-time lattices are now often used for 
this purpose in QCD. Essential results have been thus obtained. Nevertheless, these calculations 
are very laborious and have a low accuracy; in addition, it is generally difficult to estimate the 
calculation error theoretically. Therefore, it is interesting to seek other possible approaches 
to this problem. Even limited progress in this direction would allow comparing the results 
obtained by different methods. 

Long before the advent of QCD, studying the pion-nucleon interaction had been attempted 
by solving the Schrodinger equation in the Lorentz coordinate system in the framework of the 
quantum theory of pion and nucleon fields. The states were described by the method previ- 
ously found by V. A. Fock |]T[ in terms of vectors in the space that now bears his name. In 
this case, the mathematical vacuum of the Fock space coincided with the free theory vacuum. 
This approach, now known as the Tamm-Dancoff method [0, |3|, did not lead to success. The 
primary reason was the complexity of the physical vacuum state, which did not coincide with 
the mathematical vacuum. Without a description of the physical vacuum, it was impossible 
to investigate any other states. If ultraviolet (UV) and infrared (IR) cutoffs are introduced to 
make the number of degrees of freedom finite, it would be possible, in principle, to represent 
the physical vacuum as a vector of the Fock space with the free theory vacuum. However, such 
a representation proves extremely complicated because it is necessary to provide the trans- 
lation invariance of the physical vacuum and to satisfy the "cluster decomposition property" 
for vacuum expectation values. For these reasons, the Schrodinger equation in the Lorentz 
coordinates, where the evolution occurs in the conventional time, is hardly applicable to any 
calculations in the quantum field theory with a large coupling constant. 

As early as 1949, Dirac suggested a method for avoiding the difficulties related to the 
description of the physical vacuum state ||. He proposed using the light-front coordinates 
x ± = (x° x 1 , x 2 , where the Lorentz coordinates, and treating x + 

as time. In this approach, a theory is canonically quantized on the surface x + = const, and the 
generator P + of the shift along the x + axis plays the role of the Hamiltonian H . In addition, 
the generator of the shift along the x~ axis, i.e., the momentum operator P_, does not shift 
the surface x + = const, where the quantization is performed, and is kinematic, according to 
the Dirac terminology (in contrast to the dynamic generator P+). Therefore, the structure of 
the operator P_ in a theory with interaction is the same as in a theory without interaction, 
i.e., the operator P_ is always quadratic in the creation and annihilation operators a n + (p_,p_i_) 
and a n (p_,p±) and, as a rule, has the form (after normal ordering) 



where p± = (pi,P2) and the index n enumerates the species of the creation and annihilation 
operators. According to the spectral condition, the operator P_ is positive definite; therefore, 
the integration over p_ in the given formula is performed only from to oo. The operator P_ 
vanishes on the physical vacuum Q, i.e., P_|0) = 0, whence it follows that a n (p-,p±)\Q) = 0. 

For this reason, the physical vacuum \Q) can be taken as the mathematical vacuum of the 
Fock space generated by the operators a n + . No question of describing the physical vacuum 
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structure arises. The spectrum of the bound states can be found by solving the Schrodinger 
equation 



P + |*)=*/ + |tt) 

under the conditions 

p_\q?) =p'_|tf), Pj.|*) = 0, 
where p'_ are numbers and the mass squared is defined by 

m 2 = 2p' + p'_. 

Here, is a vector in the just mentioned Fock space. It is easy to satisfy the conditions 
P-\ty) = p'-\^) and P±\^) = 0, where p'_ is chosen arbitrarily. The problem is to solve 
the Schrodinger equation. Such an approach is usually called the "light-front Hamiltonian" 
approach. Obviously, it can also be used to calculate the scattering matrix. 

The scheme described faces considerable problems due to divergences that arise and the 
lack of the explicit Lorentz invariance. But the possibility of avoiding the direct description of 
the physical vacuum structure is such a considerable advantage that this approach continues to 
attract attention. Interest in it has increased with the advent of QCD. This paper is devoted 
to developing this approach. 

We note that the light-front coordinates and the similar "infinitely large- momentum system" 
are used in quantum field theory not only in the framework of the Hamiltonian approach based 
on directly solving the Schrodinger equation. Many results have been obtained by studying 
the limiting case of fast-moving reference frames in the framework of the explicitly Lorentz- 
invariant. theory of the scattering matrix or of the Green's functions [|, ||, [7], [J. But we 
consider only the Hamiltonian approach in this paper. 

For the theory with the given Lorentz-invariant initial action, constructing the light-front 
Hamiltonian P + proves a difficult problem. The primary reason is the specific divergences at 
zero values of tin' momentum p_ of virtual quanta. In particular, the invariant volume element 
on the hyperboloid p^ = m? has the form d 2 p±dp_/p^ and contains p_ in the denominator. 
The situation becomes more complicated in gauge theories. Even in the first papers on this 
problem [|K], [Uj , it became clear that canonical quantization on the surface x + = const could be 
performed only in the gauge = or in similar gauges (for example, d-A^ = 0). The reason 
is that second-class constraints arise in the theory and solving them, in particular, requires 
inverting the covariant derivative D_ = <9_ + igA_. But in the gauge A_ = 0, the Feynman 
propagator has an extra term p_ in the denominator, which strengthens the singularities at 
p_ = (at least, in the perturbation theory). 

As a consequence, a special regularization is required, which consists in "cutting out" the 
neighborhood of p- = in the momentum space one way or another and which results in 
breaking the Lorentz invariance (until the regularization is removed). This is unavoidable in 
the liglit-front Hamiltonian approach. In principle, we can preserve the gauge invariance if, 
instead of "cutting out" the neighborhood of p_ = 0, we restrict the space-time with respect 
to the coordinate x~ (— L < x~ < L) and impose periodic boundary conditions in x~ on all 



fields ||12|| . In this case, the spectrum of the momentum P_ becomes discrete, and the "zero 
modes" of the fields A, i.e., the Fourier modes corresponding to p- = 0, are explicitly singled 
out. To preserve the gauge invariance, these zero modes must be taken into account. It was 
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shown in that secondary second-class constraints arise in the theory, from which we must 
find the above-mentioned zero modes as functions of the other modes and substitute them in 
the Hamiltonian. These constraints arc so complicated that solving them proves impossible. 
Therefore, from the very beginning, we must discard the zero modes in the Lagrangian, which 
results in breaking the gauge invariance. Therefore, the breakdown of the gauge invariance by 
a regularization is unavoidable in any case. In what follows, to regularize the theory, we "cut 
out" the neighborhood of p_ = 0. In addition, the conventional UV regularization is, of course, 
necessary. 

It follows from the above discussion that the formal canonical quantization on the hyper- 
surface x + = const can produce a Hamiltonian corresponding to a theory that is not equivalent 
to the initial Lorentz-invariant one, even in the limit of removing the regularization. As a rule, 
the equivalence can be provided only by adding nonconventional counterterms |I3|, |L|, |15| to 



the light-front Hamiltonian. 

In the last few years, a number of papers have appeared in which the approximate regularized 



light-front QCD Hamiltonian was directly constructed [|16|, [Xt|] . Using the renormalization 
group theory, the form of the Hamiltonian was adjusted based on the requirement that the 
result be weakly dependent on the UV cutoff, with the relation to the conventional Lorentz- 
invariant theory not being traced in detail. Methods for simplifying the Hamiltonian and 
solving the Schrodinger equation were also proposed. The technique described in these papers 
is of considerable interest and continues to develop. But it remains unclear to what extent the 
light-front theory thus obtained corresponds to the conventional Lorentz-invariant QCD. 

In this connection, we meet the problem of constructing the light-front Hamiltonian such 
that it produces a theory equivalent to the conventional Lorentz-invariant one in the limit of 
removing the regularization. As a necessary condition, we should first provide this equivalence 
in the framework of the perturbation theory. Then, in particular, we can use the technique 



described in [16|, [L7j to simplify the obtained Hamiltonian and the subsequent nonperturba- 
tive calculations based on the Schrodinger equation. In doing so, we see the relation to the 
conventional Lorentz-invariant theory. 

The problem of constructing the counterterms for the light-front Hamiltonian, which pro- 
vide the equivalence of this approach to the Lorentz-invariant one, was investigated in in 
the framework of the perturbation theory in the coupling constant. In that paper, the authors 
proposed a method for comparing two perturbation series for the Green's functions constructed 
using the light-front Hamiltonian for one and the conventional Lorentz-invariant approach for 
the other. In addition to the UV regularization, they regularized the singularities at p_ = 
using the cutoff > e > 0, i.e., eliminating the Fourier components with \p_\ < e for every 
field in the theory. It was revealed that for the required equivalence for the nongauge field the- 
ories (in particular, for the Yukawa model), it was only necessary to add a few counterterms to 
the canonical light-front Hamiltonian. But for gauge theories (both Abelian and non-Abelian) 
under the given regularization, it proved necessary to introduce an infinite number of coun- 
terterms in the Hamiltonian. This situation is closely related to the gauge condition A_ = 0, 



which is unavoidable in the canonical light-front quantization. It is known |T8|, [[9] that to 
correctly construct the perturbation theory in this gauge, it is necessary to use the gauge field 
propagator in the form proposed by Mandelstam and Leibbrandt pX pi 



-i5 ab ( k^riy + k„n f± n k + \ -i5 ab ( fc M re„ + k v n^ {kn* 



P + iOV 2k + k„+i0 n+l P + V 2(kn*)(kn)/(nn*) + iO (nn* 
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where fi, v — +, — , 1, 2; k± = {k$ ± k 3 ) / \/2; 2 = 0, n* + 1 2 = 0. The additional pole in k_ 
in this expression is cut out with the regularization \k-\ > e. The distortion arising in this 
case does not disappear in the limit e —>■ 0. An infinite number of counterterms is required to 
compensate this distortion. 

The aim of this paper is to overcome this difficulty and obtain the required Hamiltonian 
with a finite number of counterterms. In such a case, we must change the regularization. For 
this purpose, we propose shifting the pole with respect to k_ in expression ([]]) from the point 
k_ = by changing the Lagrangian such that a small regularizing mass parameter fi 2 is added 
to the quantity 2k + k_. In this case, the distortions caused by cutting out the interval \k_\ < e 
are not so large as in the preceding case, and a finite number of counterterms is sufficient to 
obtain the correct Hamiltonian in the limit of removing the regularization (e — > and then 
fi — > together with A — > oo, where A is a UV regularization parameter). We choose the 
UV regularization such that, after removing the IR regularization (i.e., after setting e = and 
fj, = 0) at the intermediate stages, we obtain the Lorentz-invariant Lagrangian regularized in the 
UV region. This increases the number of "ghosts," but the number of necessary counterterms 
would otherwise increase sharply. We note that in the gauge A_ = 0, the one-particle irreducible 
vertex parts with the upper index "— " do not contribute to the Green's functions. In what 
follows, we prove the coincidence of every order of the perturbation theory in the limit of 
removing the regularization only for the Green's functions and not for the vertex parts with 
the index " — " . This is sufficient because the masses of the bound states are determined by the 
poles of the Green's functions. 

Under the given regularization, the Hamiltonian acts on the space with an indefinite metric, 
which prevents using the conventional variational principle to solve the Schrodinger equation 
under the conditions of preserving the regularization. Nevertheless, there exist different varia- 
tional methods that allow solving this equation. 

The chosen regularization breaks the local gauge invariance but preserves the global SU(3) 
invariance. For this reason, the number of necessary counterterms, being finite, is essentially 
larger than the conventional one. We show that there is a way to choose the counterterms such 
that the obtained light-front theory is perturbatively equivalent to the initial Lorentz-invariant 
one in the limit of removing the regularization. 

We achieve this goal by starting with the Lagrangian 



L = L + L! 



(2) 




(3) 





(4) 
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where 



fa _ q ao, _ a 40 J_ _ j V^- 2 ) 3 — 0; (K\ 

- O u A jtlx , A 2 - j 1/A 2 + 1//i2; ^ = 1; 10 ) 



v = i, X> = o, E^ M < = °> E^ M ' = o, (6) 

1=0 1=0 1=0 

A l= E ^ V> = X>- (7) 

.7=0,1 1=0 



Here 



&H = 2<9+<9_, d\ = d\ + dl N al3 



n a n */3 



7 M are the Dirac matrices, A a are the matrices of the fundamental representation of the gauge 
SU(3) group, and 

TrA a = 0, X a \ b = if abc X c + d abc X c + -5 ab , Tr(A a A 6 ) = 25 ab . (9) 

3 

We assume that the fields A* „, are restricted by the condition 

A}_ = 0. (10) 

In addition, we introduce the cutoff in the momenta k- and k± 

e<\k_\<V, v 2 <k\<V 2 (11) 

for all fields from the very beginning. This cutoff as well as condition ( |T0"D excludes certain 
degrees of freedom directly in the Lagrangian; such a procedure does not lead to new constraints 
in the canonical formalism, because no variation with respect to the excluded degrees of freedom 
is performed. 

The quantities e, A, /i, v, V and M1-M3 are regularization parameters. The coefficients 
Cj are renormalization constants, i.e., they are functions of the regularization parameters and 
are expansions in the coupling constant g. These expansions begin with g (with the coefficient 
one) for C3 and Cg and with g 2 and higher for the others. Expression (^) contains the con- 
ventional QCD interaction and the additional terms necessary for the renormalizability under 
the assumption that the Lorentz invariance and the global, but not local, gauge invariance are 
preserved. 

In Sec. [3], we prove that the light-front Hamiltonian corresponds to the given Lagrangian. 
In the framework of the perturbation theory with respect to the coupling constant g and with 
a certain dependence of the renormalization constants on the regularization parameters, this 
Hamiltonian produces the Green's functions of the fields Aq , ipo, and ipo coinciding with 
the Green's functions of the conventional QCD (renormalized in the Lorentz coordinates us- 
ing dimensional regularization) in every order with respect to g in the limit of removing the 
regularization (according to the special prescription). 
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2. The light-front Hamiltonian 

We develop the canonical light-front formalism for the theory with Lagrangian (§)-(§]) de- 
fined on the fields satisfying the condition A a -_ = and subject to cutting out the momenta 
according to formula (pi]). We use the following representation for the bispinors ipi and the 
matrices 7 M : 

Here and in what follows, the indices k, I ranges over 1, 2. 

Let us rewrite the expression for the L in a form more convenient for the canonical light-front 
formalism development using instead of the A°j + new variables 

y$ = - d k A a hk . (13) 

For example, for the Lj we get: 

L, = c' 01 Ald-d + A a k - coAldlAl + c^V + c n (^ a + d k A a k ) ^ + 
+ci 2 (</ + d k A a k ) (<p a + diAf) - c 2 A a k A a k + c 3 f abc A*A b u d»A c >» + c 31 f abc A^A b u N a » d a A c > v + 

+A«A b u A°A d s (c 4 f abe f cde g^g'' s + 5 ab 5 cd (c 5 g^g" s + c 6 g^g^)) + 
+c' 71 iV2tjj+d + tl) + + c' 72 i\^2iptd^ip^ + c 7 (ip^id±ij}^ + h.cj — c 8 (iip^ip + + /i.c) + 

+c' gi V2A a + tPX^ + + c 9 (V»ii ± ^_ + h.c.) , (14) 

where 

c 01 = ^Co + Coi, C 71 = C7 + C71, c' 72 = C 7 + C72, Cg 1 = C9 + C91. (15) 

We transform the part of the Lagrangian that contains the derivatives d + of the fermion fields 
ip, presenting it as a diagonalized bilinear form: 

1=0 Vl 1=0 1=0 Wl 

where 

V^ = £tW+, ( 17 ) 

Z'=0 

1/w/ are the eigenvalues of the matrix .B^, and t/jj/ is the matrix of the transformation diago- 
nalizing Bui. 

Consider the part of the Dirac equation that is a result of the variation with respect to ip^_: 
^ (V2d_Vz- + d ± ij l+ - M^ l+ ) + c' 72 V2d_ij„ + c 7 d ± ij + - c 8 ^+ - zc 9 A ± tP + = 0. (18) 
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This equation does not contain derivatives with respect to x + . It is a constraint that can be 
used to express ipi- in terms of other variables. To do this we sum the equation ([18]) over index 
I with the weight vi. Using the equations (0) we obtain 



3 

E 

l'=0 



\<h -M v )^v+. (19) 
By the substitution of this expression to the equation ( |i8|) one can obtain easily from the latter 

(20) 



V2 



where 



Yi = - (d x - M t ) fa+ + 



vi 



-72 



(d± - Mi>) tj) V+ - c 7 d±tfj + + c 8 ^ + + ic 9 A±tjj + 



l'=0 



(21) 



Let us remark that such a constraint can be resolved in this way not only in v4_ = gauge 
but in any gauge, if Pauli-Villars fermions are present (if no Pauli-Villars fermions we need to 
invert the operator D_ = cL — igA_, and therefore to use the A_ = gauge). 

In turn, the ip + is expressed in terms of the ip' l+ by (|7|) and (fP7|). After substituting the 
expression ( pOl) into Lagrangian the latter depends only on the variables A - and 

We eliminate the derivatives d^A°j k in the Lagrangian by integrating by parts and find the 
momenta conjugate to the A a - k (by corresponding variation of the L at fixed </?"): 



n? 



d(d + Al 



:-i) j 



Ad 2 _d + A a hk 



dld-Al k + A*(l-^ 



- c' 01 d.A a k . (22) 



If we similarly define the momenta conjugate to (p°, we obtain the second-class constraint. 
Using the Legendre transformation with respect to the variables d + A a - k and n" fc , which does 
not affect the variable <9+y?" , we pass to the first-order Lagrangian. This Lagrangian depends on 
the derivatives d + ip' h d+ip?, and d + A^ k only linearly (the dependence on d+ipl and d + ip°j is linear 
from the beginning). Using the Fourier transformation- type formulas, we then pass to the new 
variables, playing the role of the creation and annihilation operators in order that the following 
conditions be satisfied. First, the part of the Lagrangian containing derivatives with respect 
to x + must have a standard canonical form. This automatically solves the above-mentioned 
constraint for the variable (fj. Second, the positive-definite free-part must be diagonal in the 
corresponding Fock space. The following changes of variables meet these conditions: 



(2tt) 3 / 2 



J d P- J d 2 P± 



a°- rk (p)e rpx + h.c. 



e<p-<V v 2 <p 2 .<V 2 



r=0,l 



(2tt) 3 / 2 



dp- 



e<p-<V 



v 2 <p 2 <V' 



d 2 P± E ( 

r=0,l 



:(p)e" 



■ipx 



h.c. 



(2tt) 3 / 2 

£<P-<V 



d 2 p± 



v 2 <p 2 <V 2 



cij(p)e vpx — h.c. 



2 -l/4 



l+,s 



(2tt) 3 / 2 



dp_ 



d 2 p± (bUpy~ ipX + dl + (p)e 



ipx 



(23) 



£<P-<V v 2< p 2<V 2 
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where Uj = p. 



1 — jjp = p- {^jp- — I J , 0{p) = 0(j>-,px), px = p-X + pkX k , the index s 

enumerates the spinor components (s = 1, 2), and i is the index of the fundamental representa- 
tion of the color group. All creation and annihilation operators are defined for p->£. We also 
assume that v > fi in these formulas. The resulting Lagrangian has the form 



» / dP- J d 2 p ± Lj2(-iyaf(p)d + a<;(p) + 

e<P-<V v 2 <p 2 , <V 2 ^ j ' a 



+ E (-^ a jr ,k + (p)d + a% k (p) + E - {bV(p)d + bl(p) + di + (p)d + dl(p)) - H, (24) 

j,r,a,k l,i ' ) 

where H = P + is a Hamiltonian. 

Accordingly, the (anti) commutation relations have the form 



a° r ,k(P)> a h',k> + (p')} = {-l) r 5 M '8 jr 5 rr ,8 w 5{p^ - p'_)6®(p ± - p' ± ), 
[a](p),ay(p')] = -(-iy5 aa '5 jjl 5(p^-p'_)5^(p ± -p' ± ), 

{b\MX</(p')} = {d\M4/(p')} = 

= wi^'SwS^S^ -p'_)s ( - 2 \ p± - p ' ± ). 



(25) 



The negative signs in the right-hand side of the (anti) commutation relations correspond to the 
degrees of freedom carrying an indefinite metric in the space of states, and the corresponding 
operators are "ghosts". 

In contrast to the conventional canonical light-front formalism, this formalism contains no 
constraints that halve the number of the Fourier components of the fields A°- k . To preserve the 
positivity of the momentum p_ everywhere, we are forced to double the number of the creation 
and annihilation operators a°j r k + and a a - rk by introducing the index r = 0,1. The first-order 
part of the Lagrangian, which contains no derivatives with respect to x + , coincides with the 
Hamiltonian H = P + up to a sign. It has the form 



H = J dx- J d 2 x ± e 



;-iyAj 



i=o,i 



d 2 d 2 
-l) 3 (dJ-N-^ + c' 01 A k j + -jgAj >k - (l - -^) A j,h 



+ 



2 tf(l-jt)<P-\+caA k diA k + 



+4f E (9 2 ± - M 2 ) dZ%+ - -U' 72 £ (dx + Mi) dZ 1 (d x - M l ) ^+ + (26) 
v £ i=o v i y 1 i,v =o 

+ ^^^ai 1 £ (<9j_ - Mi) V>,+ + /i.c^j + (^cs^Sl 1 £ (<9± - M,) Vz+ + fr.c.^ + 
+V2c' 91 ( V a + d k A a k ) dZ 1 (V+ + (^f c^+ixcT 1 E (£l - M,) Vj+ + h.c)j - Ll I . 



Here, denotes expression ( jTjj ) with the terms with the coefficients Co, c' 01 , C7, c' 71 , c 72 , eg, 
c g , c' 91 omitted, ^+ is expressed in terms of by formulas (0) and fll7|) , and the quantities 
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A°j k , Hj k , <£j, ip' l+ and ipl^ are expressed in terms of the creation and annihilation operators by 
formulas (|23T ). 

The operator P_ has the form 

/ dp- J d 2 P± pJj2(- l ) J [- aa j + (pK(p) + 



P- I dp 

£< P -<V „2< p 2<y2 



J'=0,1 



r=0,l 



1 



+ 



1=0 Wl ) 

This operator is positive definite in the Fock space with the vacuum defined by 

o?|0) = o^|0> = ^|0> = ^|0>=0. 



(27) 



(28) 



In the framework of the perturbation theory, in the limit of removing the regularization, 
all "ghosts" are switched off in the sense that the Green's functions of this theory tend to the 
Green's functions of the correct theory, which has no "ghosts". This gives us hope that in the 
limit of removing the regularization, the unitarity condition also holds beyond the scope of the 
perturbation theory. 



3. Comparison of the light-front and Lorentz-invariant 
perturbation theories 

By the result of the perturbation theory, we mean the set of the Green's functions for the 
fields Aq , {jj , and ip constructed perturbatively in the coupling constant g. We regard the 
fields A^ , ipi and ipi for I — 1, 2, 3 as auxiliary. We show that Hamiltonian (|26|) with a certain 
A dependence of the coefficients q produces a perturbation theory coinciding in the limit of 
removing the regularization with the renormalized perturbation theory obtained from the con- 
ventional QCD Lagrangian in the gauge A- = using the Mandelstam-Leibbrandt prescription 
and dimensional regularization (the conventional perturbation theory in what follows). (See 



18] , |19f for the renormalization of the conventional perturbation theory.) The regularization is 
removed as follows: first, V — > oo, then e — > 0, then A — >• oo; Mi, M 2 , M 3 , /x and v are assumed 
to be functions of A such that Mi, M2, M3 — > 00, \x — > and v —>■ as A — > 00. The latter two 
functions must satisfy more exact restrictions: v > \x (this condition was used in constructing 
the Hamiltonian) and fiA,vA — > (these restrictions are obtained below). 

A^oo 

First, the noncovariant perturbation theory produced by the Hamiltonian can be obtained 
from the Feynman perturbation theory constructed based on the Lagrangian corresponding to 
the given Hamiltonian by resumming diagrams in every order and using the following integra- 
tion rule: in calculating the diagrams, we first integrate over k + (the momentum component 



conjugate to the light-front time x + ) and then over the other components |22, ^3fl . There- 
fore, it is sufficient to prove that the perturbation theory obtained from Lagrangian (@)-(|||) 
and supplemented by the given integration rule and the limiting transitions coincides with the 
conventional perturbation theory. 
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We represent the free part of the Lagrangian, Eq. (^), in a form convenient for the pertur- 
bation theory analysis, 

L o = I E(-i)^> (i + § - §) 0r<9 2 - &>p) a\ v + 
3 i - 

+ ]r-<M*7 M <9 M -M0^, (29) 
i=o Vl 

and we separate the part corresponding to the conventional interaction from (Q), 

L conv = g f*bc A a A b^ A c,v _ johe jcde J^^CH A <l,u _ gA a ^'f^-l(). (30) 

We assume that the remaining part of (Q) (where all coefficients are the expansions in g starting 
from the terms of order g 2 ) consists of the renormalization counterterms, i.e., it eliminates the 
divergences arising in the perturbation theory as A — > oo. The notation used and the additional 
conditions adopted are given by formulas (|5|)-(|TTD and in the text following them. 
The propagators of the fields A°j and ipi in the momentum space are 



Aab ~it ab ( n k p n v + k v n p \ (-1)' f . 



7? + ~ z0 

J 



At = Wl k*-M?+iO > (32) 

where n + = 1 and n_, n± = 0. 

Because the fields A°j and ipi always enter interaction (f|) in terms of the sums A^ and ip 
the sums of the propagators 



** = Y,*&» and A ^E A f 



always enter the diagrams as in the Pauli-Villars regularization. After all diagrams are presented 
in terms of the summary propagators, we can take the limit V —> oo for an arbitrary diagram 
(i.e., remove the restrictions \k_\ < V and k\ < V 2 ) because, in view of conditions (0), the 
propagator A^ decreases sufficiently fast (the sufficiently fast decrease of propagators fl3~T|) is 
provided by the finiteness of A) and the integrals converge. 
The summary propagator A a p b u is 



^ - k 2 + to \k 2 -p + io 9pu k 2 -fi 2 + i0 Zk+ 1 K ' 



where 



A 2 = ^7Fr4» R= T. 5 L r . „ 3V (34) 



A 2 + /^ 2 ' (l-||-i0)(l + ^ 
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After removing the cutoff (A — > oo), we have jl — > 0, and the propagator takes the conventional 
form. In terms of the propagators A^ and A^ and the vertices from fl30f ), the set of Feynman 
diagrams is the same as in the conventional perturbation theory. 

It is easy to see that as long as A is finite, there are no UV divergences in the perturbation 
theory constructed using the Lagrangian under consideration. Using this fact, as well as the 
condition k\ > v 2 (see (0)), we can apply the formalism presented in Jl3| to the perturbation 



theory and show that for the majority of diagrams, after the limit e — > is taken, the result of 
their light-front calculation (where we first integrate over k + and then over the other compo- 
nents according to the rules providing the correspondence with the noncovariant perturbation 
theory as explained above) coincides with the result of calculating the same diagrams in the 
Lorentz coordinates. The possible discrepancy that arises in some cases can be compensated 
by redefining the coefficient C2 in the Lagrangian. This is shown in Appendix 1. Therefore, 
it is sufficient to prove that in the limit A — > oo, the perturbation theory with the summary 
propagators A^J, and A^, with interaction (^), and with the restrictions e < \k-\ < V and 
k\ < V 2 removed coincides with the conventional perturbation theory. We emphasize that we 
can now perform all calculations in the Lorentz coordinates; therefore, we can make the Wick 
rotation and pass to calculating the diagrams in the Euclidean space (the location of the poles 
of propagators ( |3"3"| ) allows this). 

Propagator (^) differs from the propagator of the conventional perturbation theory, first, 
in the factor R providing the UV regularization, second, in the quantity jl, which vanishes as 
A —>■ oo and, third, in the condition k\ > v 2 , where v — > as A — > oo. 

We now analyze the behavior of an arbitrary Feynman diagram as — > and v — > (for 
finite A). In this case, essential IR divergences (essential in the sense that they arise for any 
values of external momenta and not just for special values) can occur. If such a divergence 
does not appear, then in investigating the limit A — ► for an arbitrary diagram, we can at 
once set jl — and v — in its integrand. In this case, the error in the integrand contains 
the factor jl 2 . The UV divergence of the initial diagram is not worse than quadratic; therefore, 
after separating the factor jl 2 , where ji 2 ~ /i 2 (l + k^/A 2 ), the divergence is not worse than 
logarithmic, and the condition for an error decrease is /z 2 lnA — > 0. This consideration does 

not take an increase in the IR divergence after separating the factor jl 2 into account. Its power 
increases by two. Because there was no divergence before, this power becomes not, greater than 
one, i.e., integrating the IR divergence gives (in view of the factor jl 2 ) the order jl. Integrating 
over the remaining variables produces an UV divergence not worse than linear. Consequently, 
the condition for the error decrease is fiA — > 0. Similar considerations give the condition 

A^oo 

vA — ► 0. 

We analyze the occurrence of essential IR divergences as fi, v — > in Appendix 2. It turns 
out that such divergences arise only in one case - for the diagram terms which contain the 
following factor 

A +p G pu A va , (35) 

where G pu (k) is an arbitrary one-particle irreducible two-point subdiagram, and A are propa- 
gators. In which connection the divergence takes place only when index /i in formula fl35|) takes 
the values 1, 2, and the divergence the whole of diagram at that is logarithmic: (ln/x)^, where 
is not larger then the number of factors of form (|3"5j). 

It is clear why such a divergence does not produce any problems in the conventional pertur- 
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bation theory, when G pu {k) = G p ^ irn {k) is calculated gauge invariantly with the aid of dimen- 
sional regularization. In the expression (|3^) at p = " — " the first propagator turns to zero, 
and, besides, g ++ = g + ± = 0, hence, one should consider in the propagator only the item con- 
taining the sum (k + n p + k p n + ). The first term of this sum does not give a divergence because 
n\ = n 2 = 0, and the second term gives the factor k p G pv [k\ which is equal to zero because 
of the Ward identities (analogous to ones adduced in work [25R ) which are the consequences of 
exact maintenance of gauge invariance. From the given reasoning it is clear that breaking the 
gauge invariance without the simultaneous regularization of the essential IR divergences makes 
the perturbation theory senseless. In our case, this regularization is provided by introducing 
the quantity p. 

By induction on the loop number, we prove that with a certain choice of the coefficients of 
the counterterms in Lagrangian (H) in every order, the value of every Feynman diagram tends to 
its conventional value calculated using dimensional regularization as A — ► oo. It is clear that in 
the one-loop order, there are no subdiagrams; therefore, there are no essential IR divergences. 
Then on considering a diagram containing factor fl35|) , for the lower order G pu (k) subdiagram, 
we have 

G pv {k) - G% m {k) = O (T) (36) 



A 

at p, v ^ " — " (just this estimate is obtained in Appendix 4). Because the value G pu {k) at 
p = " — " or v = " — " does not give a contribution to the expression fl3~5p then taking into 
account ( |31S| ) we can maintain that accurate to O (jA the expression (^) coincides with 



A+ p G dim A ua , (37) 

where, as it was already said, the divergence is absent. 

Therefore if (\np) N /A — > then in the limit A — > oo the diagrams containing the factor 

A^oo 

of form (^) will not differ from their values calculated in the conventional perturbation theory 
under dimensional regularization. More exactly a possible difference is due to the diagram 
divergence, but it is polynomial and is compensated by the counterterms of the same form 
that arise under the renormalization. Therefore, it is now sufficient to prove that in the limit 
A — > oo, the Euclidean perturbation theory with propagator (|33"D , where we set p = 0, with 
the propagator A^, with interaction (^), and with restrictions (|TT|) removed coincides with the 
conventional perturbation theory. 

Propagator ( |33"D with p = after the transition to the Euclidean space can be written down 

as 

- (kn) ) TT§ • (38) 

This expression is Lorentz-invariant if we assume that the vector n v (complex in the Euclidean 
space, such that in — n 3 = 0, — 0) to be properly transformed under the Lorentz transfor- 
mation. It is interesting that there is no distinguished vector other than n u in the Euclidean 
space, whereas in the pseudo-Euclidean space, the Man deist am-Leibbrandt prescription distin- 
guishes the additional surface k + = 0. However the vector n u is complex, and there is new fixed 
vector, namely, the complex conjugated to n u . It is seen that it coincides (up to a factor) with 
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Euclidean continuation of the vector n*, which is defined below the equation ([I]) and picks out 
the surface k + = in pseudoeuclidean space. 

Interaction Lagrangian fl3"0|) is also Lorentz-invariant; therefore, the counterterms that must 
be added to the Lagrangian under the renormalization in every order of the perturbation theory 
are Lorentz-invariant. It is shown in the Appendix 4 that the values which it is necessary to add 
to the diagrams in order to in the limit A — > oo make them finite and coincident with their values 
calculated using dimensional regularization are polynomials with respect to external momenta 
with the coefficients containing factors N a p (see the definition of this value in formula @). We 
should take into account that in the counterterms the vector n u cannot be contracted with the 
field A u , because we consider the Feynman diagrams for the Green's functions, whose external 
lines cannot carry the upper index " — " as it gives zero after convolution with the propagators. 
It is evident that the counterterms are globally gauge invariant because the initial Lagrangian 
has this property. We now analyze the possible form of the counterterms. 

There exist logarithmically divergent diagrams with four gluon external lines. In Ap- 
pendix 3, we show that the structure of an arbitrary diagram of this type with respect to 
the labels of the gauge group can only have the form f abe f cde or 5 ab 5 cd . Taking the Lorentz 
invariance into account, we conclude that the general form of the corresponding counterterms 
is exhausted by the terms with the coefficients C4, C5, and Cq in Lagrangian (|]) and we can 
therefore replace the addition of the counterterms by a certain choice of these coefficients. 

There also exist divergent diagrams with three gluon external lines. In general, they can 
diverge linearly; however, because of the Lorentz invariance, the divergent part must contain the 
factor and the divergence is really logarithmic. We show in Appendix 3 that the structure of 
an arbitrary diagram of this type with respect to the labels of the gauge group can only have the 
form f abc . For the general form of the divergence, this gives the terms with the coefficients C3 
and C31 in (f|) (note, that the existence of the latter takes account of the volume N a p appearance 
in the counterterms), and we can therefore replace the addition of counterterms by a certain 
choice of these coefficients. In addition, there exist logarithmically divergent diagrams with 
two fermion and one gluon external lines. It is evident that the general form of the divergence 
is defined by the terms with the coefficients eg and Cgi in @ and we can therefore replace the 
addition of counterterms by a certain choice of these coefficients. 

Next, there exist divergent diagrams with two gluon external lines. They diverge quadrati- 
cally, and the general form of the divergence is given by the term C2 in (Q). But after subtracting 
the quadratic divergence, a linear divergence can remain. Because of the Lorentz invariance, 
it is really absent, and only the logarithmic divergence exists. It is evident that the general 
form of this divergence is given by terms with the coefficients Co, Cqi, ci, en, cu in (|) and 
we can therefore replace the addition of counterterms by a certain choice of these coefficients 
and coefficient C2. There also exist divergent diagrams with two fermion external lines. They 
diverge linearly, and the general form of the divergence is given by the term with the coefficient 
Cs in (|j). But after subtracting the linear divergence, a logarithmic divergence can remain. It 
is evident that the general form of this divergence is given by the term with the coefficients C7, 
C71, C72 in (||) and we can therefore replace the addition of counterterms by a choice of these 
coefficients and coefficient c 8 . 

We conclude that the perturbation theory with propagator (PB|), with the propagator 
and with interaction (f|) (we must compare such a perturbation theory with the conventional 
one) is renormalizable by renormalizing the coefficients q. It is clear that by properly adjusting 
the additions to the quantities q in every order of the perturbation theory, i.e., by manipulating 
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the finite renormalizations, in the limit A — > oo, we can achieve the coincidence of the value 
of every Feynman diagram with its conventional value (the corresponding scheme is briefly 
presented in Appendix 4). This is just what we wanted to prove. 
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INTAS within the framework of the research program of International Center of Fundamental 
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Appendix 1 

We compare the results of calculating an arbitrary Feynman diagram in the light-front 
coordinates (with the limit e — > subsequently taken) and in the Lorentz coordinates. Every 
diagram is constructed from summary propagators (|33|) and as well as from the vertices 
entering interaction Lagrangian (f|) with the conditions \k_\ > e and k\ > v 2 but without the 
conditions |/c_| < V and k\ < V 2 . 

We use the formalism presented in [ 13|| . Under the condition k\ > v 2 , which is equivalent to 



the presence of the nonzero mass in two dimensions, the form of propagator ( |33"D is admissible 
for this formalism. It is easy to see that for all diagrams, the index u n of the UV divergence with 
respect to k + and is negative; therefore, for our theory, there are no special cases described 
|i~3[| . The numerators of all integrands of the Feynman diagrams are polynomials; therefore, 



in 



we have r > and rj > (for the notation, see |13|]). The basic formula is 



w+-n + 7j), (Al.l) 

where the minimum is taken over all subdiagrams, and the required difference between the 
light-front and Lorentz-invariant calculations is of the order e a '. 

An external gluon line carrying the label "+" contributes +1 to the value of uo_ — uo + . A 
pair of external fermion lines that are connected with the continuous fermion line contributes 
— 1 (if the diagram is proportional to 7 + with respect to the labels of this pair) or +1 (if it, is 
proportional to 7~) or (in the other cases) to the value of oj_ — uj + . We can see from formula 
( ]33| ) that without considering the factors from the vertices, the summary gluon propagator of 
the Il-line contributes —2 (if the propagator carries the labels "_L_L") or —3 (if the propagator 
carries the labels "_L+" or "++") to the value of fi. The factors from the vertices carrying 
the index "+" contribute +1 to the value of /x. From formula (|3~2"D with conditions (§) taken 
into account, we see that the summary fermion propagator of the Il-line contributes —2 to the 
value of \i. This number increases if we decrease the number of the additional Pauli-Villars 
fermion fields. Analyzing this information, we obtain the general form of the diagrams with 
o < (in fact, o = 0). It is presented in Fig. 1. The conditions for these diagrams to be 
trivially dependent on external momenta are fulfilled (see |i~3l); therefore, for these diagrams, 
the required difference is 

Cg AB , (A1.2) 

where C is a constant in the limit e —>■ 0, i.e., the dependence on external momenta is absent. 
There is no possible logarithmic correction \n(e/u), because, in view of the Lorentz invariance 
in the space of k + and k_, which holds for the Lagrangian, the quantity u must behave like 
the "—"-component of a vector, but there are no such expressions. There could only be the 
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Fig. 1. The general form of the QCD diagrams for which there exists a difference between the 
light-front and the Lorentz-invariant calculations; p is an external momentum; the symbols " +" 
and "_L" at the lines denote the corresponding Lorentz index of the propagator. 



component of an external momentum, but, we already know that there is no dependence on it. 
Therefore, all the difference is compensated by adding the term of the form 

WK> ( A1 -3) 

to the Lagrangian, i.e., by redefining the coefficient C2 in formula @. 

Appendix 2 

We analyze the possibility of the occurrence of the essential IR divergences (i.e., those 
occurring at any value of the external momenta) in the Feynman diagrams in the Euclidean 
space. We consider two cases of the divergences: when integrating over only the longitudinal 
momenta (they can appear because of (the term proportional to k + / k 2 in the gluon propagator) 
and when integrating over all momenta (the factors 1/k 2 in the propagators also contribute to 
these divergences). 

First case. We study whether a divergence exists if a part (or all) of the longitudinal loop 
momenta tend to some finite values. We assume that the external momenta do not take the 
special values (where the sum of a part of them is equal to zero). In addition, we assume that 
the transverse loop momenta are such that for all of the propagator momenta, we have Q± 7^ 
(if this condition is violated, we must take the extra contribution from d 2 Q± into account; this 
is the second case). Then the divergence can arise only if for some gluon line, we have Q n = 0. 
The factor in the integrand producing the possible divergence has the form 

Q^n u + QyUp _ /^ (Q^n u + Q v n^){iQ + Q 3 ) /Aon 

Q a n° ~ " V2 qTTqI (A2 - 1} 

and is a pole of an order not higher than one. In the loop momentum space, we consider a 
point where Q n = for a certain set of lines. We look for the power a of the IR divergence. 
Every line of this set contributes —1 (if it carries the labels +_L, see (A2.1)) or (if it carries 
the labels _L_L or ++). We exclude all lines with the labels _L_L and ++; then every line of 
the set contributes — 1. The differentials d 2 q n of the loop momenta (the integration volume 
elements) give a positive contribution to a. We must consider only those loop momenta whose 
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change (other momenta being fixed) violates the condition Q n — for the lines of the set. The 
number of such loop momenta is the number of lines whose momenta can be taken arbitrarily, 
i.e., the number of independent lines. The total positive contribution to a is equal to twice this 
number. We then find other contributions to the IR divergence. 

We break all lines of the set. The diagram splits into n + 1 connected parts (n = 0, 1, . . .). 
All momenta external with respect to the whole diagram enter one part (the external momenta 
would otherwise take the special values). We call this part separated and the other parts 
nonseparated. 

A nonseparated part is a subdiagram; if it has the external Lorentz labels, then it must be 
proportional to its external line momenta carrying these labels. But the factor of proportion- 
ality cannot contribute to a, because of the Lorentz invariance, the invariance with respect to 
multiplying by a complex number, and the condition Q± ^ 0. Every external line carrying 
the label k (k = 1,2) gives the factor k k or g^ k \ the line carrying the label "+" gives the 
factor k + , where A; is a linear combination of the external momenta of the subdiagram, i.e., the 
momenta of the set lines. Therefore, every external line carrying the label "+" contributes +1 
to a, and one carrying the label k contributes zero. Let m be a summary contribution to a 
obtained in such a way from the nonseparated parts. 

Let s be the number of lines of the set and r be the number of lines of the set external with 
respect to the separated part. Every line carries the label "+" at one end; therefore, we have 

m> s-r. (A2.2) 

The number of the independent lines in the set is equal to s — n (each of the n+1 parts gives 
the ^-function, and one (^-function is common to the whole diagram). By the definition of a, 
we have 

a = 2(s-n)+m- s. (A2.3) 

Using (A2.2), we obtain 

a>2(s-n)- r. (A2.4) 

Every nonseparated part has a minimum of two external lines. This means that all parts 
(separated and nonseparated) together have a minimum of 2n + r external lines of the set. All 
these lines are pairwise connected with each other. Consequently, we have 

1 r 
s > -(2n + r) = n + -, (A2.5) 

whence we obtain a > 0, i.e., the divergence cannot be worse than logarithmic. 

We find the general case producing this divergence. In this case, all the above-cited in- 
equalities must reduce to equalities, i.e., we have m = s — r, and all the nonseparated parts are 
two-point diagrams. The general form of such a diagram is shown in Fig. 2. In this diagram 
every chain giving essential IR divergences contains the factor 

A +p GrA va , (A2.6) 

where G pv is the nearest to the separated part nonseparated one which is one-particle irreducible 
subdiagram and A are the propagators of the lines entering to the G pv . The divergence takes 
place only when the index p in (A2.6) takes values 1, 2. 
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Fig. 2. The general form of the diagram producing the essential IR divergence in QCD; the 
symbols "+" and "_L" by the lines mean that the propagator has the corresponding Lorentz 
index. 



Second case. We now assume that a part or all of the four- dimensional momenta lend 
to certain finite values. The consideration is similar. The gluon propagator now gives a pole 
of the second order, whereas the fermion propagator gives no poles, because of the presence 
of the mass. Notice, that at the investigation of IR divergences one can replace all fermion 
lines by 1/M which makes impossible getting 1/M to the numerator. Besides in this case 
one can replace R by 1 and fx 2 by fi 2 in the propagator of gluon which excludes A from the 
consideration. That is why from dimensional considerations, it follows that every nonseparated 
part contributes no less than four minus the number of its external lines to a. This is also 
true for the divergent diagrams with the stipulation that for the two-point diagram, we must 
subtract the part independent of the momentum and proportional to g^ u together with the 
divergent part, i.e., we must normalize the gluon mass to zero in every order. The latter does 
not violate the Ward identities in the limit of removing the regularization. Therefore, we now 
have 

m > An — (2s — r), 
a = 4(s-n)+m-2s, (A2.7) 

whence we obtain 

a > r > 0, (A2.8) 

i.e., there is no divergence. 

Appendix 3 

In the framework of the perturbation theory constructed based on Lagrangian we 
consider an arbitrary diagram with three gluon external lines and investigate its structure with 
respect to the labels of the global gauge group. This structure has the form A abc . The gluon 
propagator is proportional to 5 ab , the three-point gluon vertex is proportional to f abc , and the 
four-point vertex is proportional to f abe f cde or S ab S cd . In addition, there can be fermion loops. 
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Every fermion loop gives the factor 

ReTr(A a ...A c ), (A3.1) 

if the number of vertices in the loop is even or 

TmTr(A a ...A c ), (A3.2) 

if this number is odd (similar to the Furry theorem in quantum electrodynamics). We consider 
only the factors related to the gauge labels. From it follows that 

KxP^d = 2 (fiasSp-y — -fiapdyS^ ■ (A3. 3) 

We can also write 

Tr (A a A fe A c ) = i2f abc + 2d abc , (A3.4) 

whence it follows that 

f abc = -~ (Tr (A a A 6 A c ) - Tr (\ c \ b \ a )) . (A3.5) 

We replace all f abc in the diagram with the right-hand side of formula (A3. 5) and then 
replace all resulting expressions of the form A a A a according to formula (A3. 3). As a result, 
only three matrices A carrying three external labels of the diagram remain. These matrices 
are connected with each other by their labels of the fundamental representation (because after 
formula (A3. 3) is used, only the Kronecker symbols remain), i.e., A abc consists of terms of the 
form 

Tr(A a A b A c ). (A3.6) 

After using formula (A3. 5), we have the factor % , where N is the number of three-point 
gluon vertices. The use of formula (A3. 3) gives no new imaginary factors. The initial expression 
consists of the real quantities 5 ab and f abc and the quantities (A3.1) and (A3. 2), of which the 
first is real and the second is imaginary; therefore, it is proportional to i N , where N' is the 
number of fermion loops with an odd number of vertices. It is easy to show that for the diagram 
with three gluon external lines, the quantities N and N' have different parities. Therefore, for 
the powers of i to be consistent, it is necessary that the sum of the expressions of form (A3. 6) 
be imaginary, whence, using (A3. 4), we conclude that 

A abc _ jabs^ ( A3 J ) 

We now similarly consider the structure of an arbitrary diagram with four gluon external 
lines. It has the form A abcd . The quantity A abcd consists of terms of the forms Tr(A a A 6 A c A d ) and 
Tr(A a A fe )Tr(A c A d ) (the latter expression is explicitly real). It is easy to show that for the diagram 
with four gluon external lines, the quantities N and N' have the same parities. Therefore, the 
sum of the indicated expressions, which compose A abcd , must be real. We can represent the 
quantity Tr(A a A 6 A c A d )) as a sum of its symmetrized part (explicitly real) and expressions of 
the form if abe Tr(\ c \ d \ e ). In view of the reality condition, the latter expression must be a sum 
of quantities of the form f abe f cde _ Direct calculation shows that the symmetrized part of the 
quantity Tr(A a A 6 A c A d ) is proportional to the symmetrized part of the quantity 5 ab 5 cd . We can 
therefore conclude that A abcd consists of terms of the forms 
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Appendix 4 



We shall find a form of the values which are necessary to add to Feynman diagrams in 
order to in the limit A — > oo make them finite and coincident with its value calculated using 
dimensional regularization. We consider a diagram in Euclidean space constructed from prop- 
agator ([53), from the propagator A^, and from the vertices entering ( p?0| ) with restrictions (JTTJ) 
removed. The other vertices entering (|4]) are involved only in reducing the subdiagrams of the 
preceding orders of the perturbation theory to the "correct" (i.e., calculated using dimensional 
regularization) value including cancellation of divergences. 

It is clear that the summary fermion propagator can be represented as 

= - k ^ R> + - M ° 9 R", (A4.1) 

k 2 - M 2 + iO k 2 — Mq + iO V ; 

where R' and R" are cutoff factors that properly decrease and R', R" — > 1. We represent the 

A^oo 

diagram as a sum such that one summand contains only one term of the numerator of every 
fermion propagator. Each of these summands can be represented as 

1 = J dkF(k,p)f A (k,p), (A4.2) 

where dk represents all volume elements, f A (k, p) is the product of all factors R' and R" entering 
the fermion propagators and of all factors R entering boson propagators (|3~3D, k denotes the 
integration momenta, and p denotes the external momenta of the diagram. 

It is evident that we can find a function F(k,p) that is polynomial in p, is Lorentz invariant, 
and has no nonintegrable IR singularities such that for F(k,p) and F(k,p), a number of the first 
terms of their asymptotic expansions at infinity with respect to k coincide and the difference 
F(k,p) — F(k,p) is integrable (see the refinement below). Therefore, we can write 

I = J dk(F(k,p) - F(k,p))f A (k,p) + J dk.F(k,p)U(k,p). (A4.3) 

Up to corrections of the order 1/A, we can neglect the factor f A (k,p) in the first integral 
(because the integral converges even in the absence of this factor). We can then assume that 
this integral is calculated using dimensional regularization and we can split it into two integrals 
(assuming that they are both renormalized by dimensional regularization). As a result, we 
obtain the expression 

/dim p-dim A /■ A 

dkF(k,p)- dkF(k,p)+ dkF(k,p)f A (k,p). (A4.4) 

Using the expansion for the function f A (k,p) in the last integral in (A4.4) with respect to p in 
the neighborhood of the origin (this expansion is well defined for any k), we can now show that 
the last integral in (A4.4) is a polynomial in p plus corrections of the order 1/A. 

We must make the following refinement of these considerations. Before subtracting the 
overall divergence of the diagram, we must verify that the diagram has no sub divergences 
with respect to a part of the integration variables. Such sub divergences can be produced by 
subdiagrams (which is taken into account by the renormalization in the lower orders) or by 
the divergence on integrating over a part of the momentum components. For the QCD in 
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the gauge A_ = 0, the latter is possible because of the improper decrease of the propagator 
in the direction k±. It is known p4| that this divergence is present if the index of the UV 
divergence with respect to a part of the components is nonnegative. From the structure of the 
propagators, we can see that in the Euclidean space, the subtraction of the overall divergence 
cannot decrease only the index u± of the UV divergence with respect to k±. This results in the 
necessity to preliminarily subtract, the subdivergence with respect to k±, the subdivergence 
being in general dependent on the projections p v n il of the external momenta in an arbitrary 
nonpolynomial way Analysis of the diagrams for the Green's functions shows that we have 
u± > for only the one- loop two-point diagrams and that u± = 0. These diagrams have only 
one external momentum. Therefore, they cannot have a nonpolynomial dependence on p u n u , 
because of the invariance with respect to multiplying the vector n v by a complex number. One 
must take into account that, for the diagrams of the Green's functions, the vector n v with the 
nonconvolute index v cannot stand in the numerator, because this gives zero on convolution 
with the propagators. If we consider the one-particle irreducible vertex parts, whose diagrams 



do not satisfy the last condition, the divergent parts can be nonpolynomial [jn|, [L9|, and the 
number of diagrams that are divergent with respect to k± can be much larger. 

It seems that if the integral of the form, described at the beginning of this Appendix, 
converges then the result of it' s calculation cannot depend on the vector n* (complex conjugated 
to n u up to a factor). But this would be true only if this integral be an analitical function of 
the complex vector n v . However the derivative of our integral with respect to complex vector 
n u can be nonconvergent due to the rising of infrared singularity (of the power of the pole in 
(nk)). Therefore the result of the calculation of the diagram, and, hence, it's divergent part 
can depend on n v and on n*. It is seen from the equation ( |38| ) for the propagator that the 
integrands we considered are invariant with respect to a multiplication of the vector n v by a 
complex number (and of the n* by complex conjugate number). This allows to conclude that up 
to corrections of the order 1/A the difference between the diagram described in the beginning 
of this appendix (more precisely, the finite sum of such diagrams of the given order) and its 
value calculated using dimensional regularization (i.e., the similar sum of the first integrals in 
(A4.4)) reduces to a polynomial in the external momenta with coefficients containing factors 
N a p (see the definition of this value in formula (P)). 
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